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Abstract. Young supernova remnants are thought to be the sites where cosmic ray acceleration occurs by the mechanism of 
diffusive shock acceleration. The maximum energy gained in this process is conventionally extimated to have a value close to, 
but distinctly below, the "knee" (~ 10'^ eV) of the cosmic-ray spectrum. Bell & Lucek (2001) have suggested that the generated 
cosmic rays simultaneously amplify the magnetic field around the supernova remnant shock to many times its pre-shock value. 
In this case the acceleration rate may be significantly increased and protons can easily reach energies up to 10'^ eV. We use a 
"simplified" box model incorporating the magnetic field amplification suggested by Bell & Lucek to investigate the resulting 
modifications of the cosmic-ray spectrum. The model predicts a spectral break at high energies, close to the "knee" region, and 
in good accordance with observations. 



1. Introduction 

Bell & Lucek (2001) and Lucek & Bell (2000) have presented 
numerical simulations suggesting that the conventional pro- 
cess of particle acceleration in shocks, generally called diffu- 
sive shock acceleration, may also result in substantial amplifi- 
cation of the highly tangled magnetic field around the shock. 
The unusual strong magnetic field inferred in Cas A (e.g. Vink 
& Laming 2003) seems to confirm this prediction. We note in 
passing that this hypothesis provides a concrete physical mech- 
anism for one of the many ideas of the late F Hoyle who spec- 
ulated (Hoyle 1960) that strong interstellar shocks might con- 
vert mechanical energy into roughly equal amounts of mag- 
netic field energy, cosmic ray energy and thermal energy. It is 
also closely related to the "plastic deformation" of the magnetic 
field discussed qualitatively by Volk and McKenzie (1981). 
Bell & Lucek (2001) suggest that this may enable strong su- 
pernova driven shocks in the interstellar medium to accelerate 
protons to energies well beyond what are conventionally held 
to be the maximum realistically attainable energies of a few 
times 10'"* eV (e.g. Lagage and Cesarsky, 1983). 

The Bell-Lucek hypothesis is of great interest because it 
is one of the few suggestions as to how the cosmic ray par- 
ticles of energies at and above the "knee", located at a few 
times lO'^ eV, could be made in the Galaxy. It is therefore im- 
portant to determine what the expected form of the source en- 
ergy spectrum would be and whether the slight break in the ob- 
served spectrum at the "knee" can be accounted for in this way. 
There are also interesting implications for gamma ray observa- 
tions of SNRs with the next generation of imaging atmospheric 
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Cherenkov telescopes such as HESS. For an initial examination 
of this problem the so-called "box" models of particle acceler- 
ation (Drury et al. 1999) offer enough accuracy and have the 
great advantage of computational simplicity. 

2. The Bell-Lucek hypothesis 

Bell & Lucek point out that in the conventional diffusive shock 
acceleration picture the energy density in the accelerated par- 
ticles at and near the shock front is extremely large, of order 
the total ram pressure of the incoming plasma, and thus much 
larger than the energy density of the magnetic field (the ratio 
of particle energy to field energy is of order the square of the 
shock Alfven Mach number). The standard treatment of reso- 
nant wave excitation, in which the perturbations of the field are 
treated as Alfven waves to lowest order and the interaction with 
the particles as a small perturbation, is thus very questionable. 
With the support of numerical simulations and simplified ana- 
lytic models they suggest that in reality the field can be highly 
distorted by the particle pressure and wound up to the point 
where approximate equipartition holds. 

If this amplified and distorted field is then used to estimate 
the particle diffusion in the shock neighbourhood, acceleration 
to substantially higher energies than in the conventional pic- 
ture is obviously possible. Detailed estimates and simple di- 
mensional analysis agree that the maximum particle rigidity 
is given, to order of magnitude, by the product of the field 
strength B, the shock radius R and the shock velocity U. Thus, 
other things being equal, the increase is directly proportional to 
the increase in the field strength, which by the above argument 
is of order the Alfven Mach number of the shock Alsh- This can 
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easily be Alsh - 10^ for a young supernova remnant so that the 
effect is potentially very significant; acceleration to rigidities 
of a few lO'^ V, rather than the 10^* V normally estimated, is 
easily possible. 

Ruskin and ZirakashviU (2003) have performed a sophisti- 
cated analysis of the combined effects of field amplification and 
wave damping on particle acceleration in supernova rermiants 
and conclude that indeed the upper cut-off to the accelerated 
spectrum can be quite strongly time-dependent. Our aim in this 
short paper is to carry out a first exploratory analysis of the ef- 
fect this has on the overall spectrum and we therefore use the 
simple approach of assuming an equipartition field although 
this is clearly a rather crude approximation. 

3. The "box" model 

To get a first estimate of the effect of such a dynamically gener- 
ated field on the acceleration we turn to the simplest treatment 
of shock acceleration, the so-called "box" model. In this the 
accelerated particles are assumed to be more or less uniformly 
distributed throughout a region extending one diffusion length 
each side of the shock, and to be accelerated upwards in mo- 
mentum space at the shock itself with an acceleration flux 
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per unit surface area where U\ and t/a are the upstream and 
downstream velocity and f{p) is the phase space density of 
the accelerated particles (assumed to have an almost isotropic 
distribution). If the diffusion length upstream is L\, and that 
downstream is L^, then 



^i(p) 



L2 « 



^2(P) 
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where ki and K2 are the upstream and downstream diffusion 
coefficients. To a first approximation we assume that both Li 
and L2 are small relative to the radius of the shock and that we 
can neglect effects of spherical geometry (in fact it is not too 
difficult to develop a spherical box model, but it unnecessarily 
complicates the argument) so that the box volume is simply 
A(Li + L2) where A is the surface area of the shock. The basic 
"box" model equation is then simply a conservation equation 
for the particles in the box; the rate at which the number in the 
box changes is given by the divergence of the acceleration flux 
in momentum space plus gains from injection and advection 
and minus advective losses to the downstream region. 

-[A(Li +L2)47rpV(p)] + A—^AQ(p) 



dp 

+ AF.ip) - AF2ip), 



(3) 



where Q(p) is a source function representing injection at the 
shock (only important at very low energies), Fi is a flux func- 
tion representing advection of pre-existing particles into the 
system from upstream (normally neglected) and F2 is the flux 
of particles advected out of the system and carried away down- 
stream. The only complication we have to consider is that the 
box is time-dependent, with flow speeds, shock area and diffu- 
sion lengths all changing. 



The escaping flux is determined simply by the advection 
across the downstream edge of the box, that is 



F2ip)=Wfip){U2-^ 



(4) 



where we have to explicitly allow for the time varying size of 
the downstream region. Substituting this expression for F2(p) 
and neglecting the advection of prexisting particles (the Fi(p) 
term) the box equation simplifies to: 



1 OA dLi df 



+ (f/l 



Q 



(5) 



4. Incorporating the Bell-Lucek effect In the "box" 
models 

Partial differential quations of this form always reduce, by the 
method of characteristics, to the integration of two ordinary 
equations, one for the characteristic curve in the {p, t) plane 



dp _ U1-U2 P 
At ~ Li +L2 3' 

and one for the variation of / along this curve 
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Apart from at the injection momentum Q = and we can write 
the above equation as 



din/ 



1 dA 

'A~dt 



1 dL, 



f/i 



df A dt L\ + Lj dt L\ + L2 
But the shock area A is a function only of time so that 
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and, although the upstream diffusion length does depend on 
both time and momentum, if we assume Bohm scaling for the 
two lengths so that 



K pv 

L oc — oc , 

U UB 

(where v is the particle velocity) we can write 
I 



Li -I- L2 dt 
where 



dLi d ln(f/i5i) 
= —ir- 



dt 



L1+L2' 

(obviously < & < I). Finally, noting that 

Ui _ 3Ui dlnp 
Li +L2 ~ U1-U2 dt ' 
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and assuming that the shock remains strong, which yields a 
fixed compression ratio, we can simplify equation ^ to 

din/ dlnA d ln(f/iBi) 3Ui d In p 

— +17- 



(14) 



df dt dt U1-U2 dt 

which integrates trivially to relate the value of / at the end of 
one of the characteristic curves, say at the point {pi, fi), to the 
value at the start, say at (fo,/7o), as follows; 



(15) 



(16) 



f(tupi) ^ iMh)] / t/i(fi)gi(fi) r Ip^ 

f(tQ,PQ) \A(to)l \Ui(to)Bi(to)j \po 

where 

3Ui 
'~ U1-U2 

is the standard exponent of the steady-state power-law spec- 
trum associated with shock acceleration. 

This rather beautiful result shows how the standard test par- 
ticle power-law is modified by a combination of effects as the 
box volume changes. As one would expect the amplitude varies 
inversely as the shock area and also decreases if the upstream 
diffusion length (at fixed energy) increases, but with an expo- 
nent between zero and one determined by the ratio of the up- 
stream diffusion length to the total width of the diffusion re- 
gion. It is very interesting that the result is not simply a vari- 
ation inversely as the box volume, which one would naively 
expect from geometrical dilution. This reflects the fundamen- 
tal asymmetry between the upstream and downstream regions, 
that upstream is empty outside the diffusion region whereas the 
entire downstream region is filled with accelerated particles. 

If we assume pure Bohm scaling the other diffential equa- 
tion is also integrable so that the problem is reduced entirely to 
quadratures (of course only within the various approximations 
we are making; but still a remarkable result). Bohm scaling im- 
plies that the mean free path is of order and proportional to the 
particle gyroradius, so that if the particle charge is e 
pv 



Li + L2 ~ a- 



(17) 
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where c is a dimensionless parameter (probably of order ten). 
Substituting in the equation of the characteristic (equation 
(|6})we get 



dp 1 

^ -(Ui-U2)UieBu 
dt a 



(18) 



and noting the relativistic identity between kinetic energy T 
momentum p and velocity v, 

dT 

we can integrate this as 



(19) 



Ti-To 



(Ui - U2)UiBidt. 



(20) 



For relativistic particles the kinetic energy and the momentum 
are essentially interchangeable with T - c^p"^ + nf-c^ - mc^ ^ 
cp. 

These two integrals (equations (I15> & ilOi ) together reduce 
the problem of calculating the final spectrum to that of deter- 
mining the initial amplitude /(fo, po) which in turn depends on 
the injection rate and its time dependence. 



5. The injection rate 

There are two main approaches to the injection rate. The sim- 
plest, which is perhaps more consistent with the test particle 
approach, is to simply parametrise it by assuming that some 
fraction of the incoming thermal particles are "injected" as non- 
thermal particles at some suitably chosen "injection momen- 
tum" which separates the thermal particle population from the 
non-thermal. In other words one writes 



Qip, t) = 77(r)ni Ux6(p - pinj(f)), 



(21) 



where n\ is the upstream thermal particle number density, rj <s: 
1 is the injection fraction, pmj is the injection momentum and as 
usual 6 is Dirac's delta distribution. It should be clear that this 
is a parametrisation rather than a true injection model, however 
it, or equivalent parametrisations, have been very widely used, 
typically with 77 taken to be a constant of order 10"^ to 10"^ 
for protons and pi„j w IQmpUi where nip is the proton mass. 
However there is no real justification for this apart from the fact 
that it seems to yield reasonable results in many cases. 

With the above parametrisation the distribution function 
just above the injection energy can be simply determined by 
equating the acceleration flux to the injection flux. 



4npl 



■(Ui-U2)f(pini)^WlU 



giving 

_ 3 Ui 

/(finj) - - —Jr 



U2 



(22) 



(23) 



The second approach adopts the idea, which can be traced 
back to the early work of Eichler, that the injection process 
is inherently extremely efficient but that various feedback pro- 
cesses operate to reduce it to the point where the accelerated 
particles carry a significant part of the energy dissipated in the 
shock. Probably the most sophisticated modern version of this 
idea is to be found in the papers by Malkov (eg Malkov 1998, 
Malkov et al 2000; see also Kang et al 2002). This, or some- 
thing similar, is in fact required for the Bell-Lucek hypothesis 
to operate because it requires the accelerated particle energy 
density to be substantial and of order the ram pressure of the 
upstream flow. For a standard spectrum close to the energy 
is almost uniformly distributed per logarithmic interval over 
the relativistic part of the spectrum. This suggests taking a ref- 
erence momentum in the mildly relativistic region, po k mc, 
and determining / by a relation of the form 



——f{po)mc 



'.l3pUi(Ui-U2) 



(24) 



where is a number which depends logarithmically on the 
upper cut-off and which for supernova remnants is probably 
somewhere between 10"' and 10"^. 

It is important to note that both injection models are mod- 
els for proton injection, the protons being the dynamically 
dominant species. Unfortunately very little is known about the 
factors controlling the injection of electrons and other minor 
species despite their importance for diagnostic tests. It is also 
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ln(p) 




ln(t) 

Fig. 1. The characteristic acceleration curves in the t, p plane 
for a Sedov blast wave and equipartition field amplification. 
The dotted line indicates the p\ oc t^^^^ asymptotic power-law 
relation between starting time and final momentum. 



very probable that the injection is nonuniform over the shock 
surface with a strong dependence on the angle between the 
mean background field and the shock normal (Volk et al, 2003). 

6. Application to the Sedov solution 

Let us now apply these ideas to the Sedov solution (also studied 
by Taylor and von Neumann) for a strong spherical explosion 
in a cold gas where the shock radius expands as /? oc fil^ and 
the shock velocity decreases as f/ oc r^l^ . On the Bell-Lucek 
hypothesis the magnetic field also scales as the shock veloc- 
ity, B oc r^l^ and thus the characteristic acceleration curves 
(eauationJ20» are given by 



T\-To oc 



f 

J to 

f 

J to 



{Ui - U2)UiBdt 



oc - fj 



4/5 



(25) 



These curves, illustrated schematically in Fig. 1, all rise ex- 
tremely steeply, representing an initial phase of rapid accelera- 
tion, turn over and then become asymptotically flat. Physically 
it is clear that, as the shock slows and the field drops, the high 
energy particles cease to be significantly accelerated and sim- 
ply diffuse further and further in front of the shock. In fact in 
reality they should probably be thought of as decoupling from 
the shock and forming part of the general interstellar popula- 
tion at this point, but within the box model they simply fill a 
steadily growing upstream region. We will return to this point 
later 

A very important aspect of the curves is that they uniquely 
relate final energies (or equivalently momenta) to starting 
times. Asymptotically the relation is a simple power-law; for 
Ti 3> To and fo f i we have simply 



pi oc Ji oc 



4/5 



to oc p 



-5/4 



(26) 



Using this we can now translate the dilution factors from equa- 
tion (I15> to additional power-law terms in the final momentum. 
Explicitly, a given final momentum maps to a starting radius 
using a Sedov expansion-law: 



R(to) cc tf oc p-'l\ 



(27) 



and thus the first term on the RHS of equation M5\ translates 
to a Pj^' factor: 

1 



A(fo) 



(28) 



The final momentum can also be mapped to a starting velocity, 
using again a Sedov expansion-law: 



uito)oct;^i'^p'^\ 



(29) 



The magnetic field, which on the Bell-Lucek hypothesis scales 
as velocity, gives an additional p J^'* factor, thus the second term 
on the RHS of equation M5\ scales as: 



f/l(fl)Bi(fi)f _3,,/2 
oc p, 

t/i(fo)Bi(fo)/ ^' 



(30) 



Furthermore we need to determine the initial amplitude of / 
from an injection model. First we use the rj parametrisation as 
discussed in section 5: 



Pmi ^ po'^ U{to) oc p^l"^, 
together with 



(31) 



(32) 



Finally assuming a strong shock which, using the Rankine- 
Hugoniot conditions for a non-relativistic fluid, yields U\IU2 - 
4 and 3Ui/(Ui - U2) = 4 and substituting the equations (128> . 
i30\ . i3H and ( I32l l into equation il5\ one obtains a scaling-law 
for the particle distribution fipi) at a fixed time t^: 



fipi) oc T^nopi^A{to)[umBmrH— 

\Po 

3/4 -1 -3»/2 -4 

Wi Pi Pi Pi ■ 



-4 



(33) 



If 77 is constant, the slope is steepened from the canonical value 
of 4 to 



3)? 

4.25 + —. 
2 



(34) 



If we use the alternative "equipartition" argument as an in- 
jection model, meaning that fo should be dynamically deter- 
mined in the mildly relativistic region, we have po ~ mc inde- 
pendent of p\ and 



f{po) oc f/f oc r ^ ocp^^^. 



(35) 



which gives the following scaling-law for the particle distribu- 
tion: 



f, ^ 3/2-1 -MI2 -A 
f(P\)^Pl Pi Pi Pi ■ 

In this approach the slope is given by 

. 3??- 1 
A+ . 



(36) 



(37) 
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Fig. 2. Example of f{ti,p\) (dotted line) from equation (15), 
where a Sedov-like shock has been to describe the velocity and 
■& - 0.9. The solid line is the standard /j"^ curve, the dotted line 
converges towards the spectral index given by equation (34). 

It is interesting that because of the strong injection at early 
times this model can even, if i? < 1/3, lead to a slight flat- 
tening of the spectrum. However, especially at high energies, it 
is unlikely that the upstream diffusion region could be so small 
and a modest steepening of the spectrum is more likely. 

These results refer of course only to the asymptotic be- 
haviour of the high energy part of the spectrum. As pi is de- 
creased there comes a point where fo is no longer small relative 
to f 1 . At this point all values of the final momentum map down 
to a small approximately constant region and the spectrum be- 
comes simply the standard p^'^ spectrum. This break occurs at 
the point to which efficient acceleration is possible at that stage 
in the remnant evolution, and decreases as the the remnant ages. 
Shock acceleration will terminate when the shock is beginning 
to weaken and the amplified field is only a few times the ambi- 
ent field, which for ambient fields of a few yuGauss and typical 
SNR parameters places the break exactly in the "knee" region. 

Figure 2 shows an example of the total spectrum for 
a Sedov-like shock, which was obtained by using equation 
(15). The injection mechanism for this example was taken ac- 
cording to the equations (22) and (23) and )? - 0.9. The spec- 
trum clearly shows the smooth transition from the standard 
spectrum to the asymptotic spectral index given by equation 
(34). 

7. Conlcusions 

We have applied the Bell-Lucek hypothesis to a Sedov-like 
shock, using a simplified Box model. We showed that such a 
model exhibits a spectral break at an energy determined by the 
current acceleration cut-off below which one observes the stan- 



dard shock acceleration spectrum, but above which a slightly 
different power-law continues to higher energies. For older 
remnants the break is expected to be in the "knee" region at 
rigidities of order 10'^ V. The key point is that the dynamical 
field amplification both increases the maximum attainable en- 
ergy and makes it a relatively strongly time-dependent quantity 
unlike the situation with no field amplification where, as is well 
known, there is only a very weak dependence of the cut-off en- 
ergy on the remnant age, at least during the Sedov phase (the 
much more sophisticated analysis by Ptuskin and Zirakashvili, 
2003, is relevant here). 

While the spectrum at the shock is the more relevant quan- 
tity from the point of view of gamma-ray tests, in the context of 
cosmic ray propagation theory what one would really like is the 
effective source spectrum of an individual supernova integrated 
over its history. This, while obviously related to the spectrum 
discussed here, is a somewhat different quantity and harder to 
evaluate. However it is clear that here also one should expect a 
broken power-law with a relatively small break in the exponent 
at a rigidity corresponding to acceleration in a mildly amplified 
Galactic field. 

Obviously further and more detailed work is needed, but 
it is very encouraging that even such a simple model can pro- 
duce spectra remarkably close to the inferred cosmic ray source 
spectrum through the "knee" region. In fact we are not aware 
of any other acceleration model that can naturally produce a 
break of the right magnitude (about 0.5 in the exponent) at the 
right position (modulo major uncertainties in interstellar prop- 
agation at these energies of course). 
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